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Abstract. We study modulation spectroscopy of the potassium D2 transitions 
at 766.7 nm. The vapour pressure, controlled by heating a commercial 
reference cell, is optimized using conventional saturated absorption spectroscopy. 
Subsequent heterodyne detection yields sub-Doppler frequency discriminants suitable 
for stabilizing lasers in experiments with cold atoms. Comparisons are made between 
spectra obtained by direct modulation of the probe beam, and those using modulation 
transfer from the pump via nonlinear mixing. Finally, suggestions are made for further 
optimization of the signals. 
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1. Introduction 

Laser cooling lias revolutionized our ability to probe and manipulate atomic gases. In 
order to perform experiments with cold atoms, lasers first need to be stabilized near 
strong cooling transitions with residual frequency fluctuations well below the natural 
linewidth (a few MHz for alkali atoms). The standard technique for obtaining sub- 
Doppler spectra in room-temperature gases is known as saturated absorption [1]. In 
this method a laser beam is divided into a relatively intense pump and weak probe, 
which counter-propagate through a gas. For simplicity, consider the case of two-level 
atoms. Since the pump and probe have equal and opposite wave- vectors, atoms with zero 
Doppler shift along the optical axis will be simultaneously resonant with both beams. 
The pump excites these atoms, reducing the ground state atomic population and causing 
reduced probe absorption in the vicinity of the resonance. The resulting transmission 
of the probe therefore shows a Doppler-broadened background with a narrow peak of 
reduced absorption known as the Lamb dip [2], [3]. 

For laser stabilization, an electronic frequency discriminant is needed with 
approximately linear slope around zero volts at the lock point [U |5]. One can lock 
directly to the side of a saturated absorption peak by subtracting an electronic offset, 
but this results in undesirable sensitivity to fluctuations and drifts in optical power, 
vapour pressure, and offset voltage. A more stable discriminant can be obtained through 
phase or frequency modulation and subsequent phase-sensitive detection. There are 
a variety of heterodyne techniques which have been used for studying sub-Doppler 
spectroscopy [6l [71 [H HI [101 [H] • These vary in the details of the optical configuration, 
the number of frequency components and how they are produced, and the relative 
strengths and detunings of the contributing beams. Beams may be generated by phase, 
frequency, or amplitude modulation, frequency shifting, or from independent sources. 
For small modulation depths, and modulation frequencies much smaller than the natural 
linewidth, phase modulation of the probe results in a signal which is approximately 
proportional to the derivative of the absorption feature [12]. For large modulation 
frequencies, spectra exhibit multiplet structures, similar to those obtained with the 
popular Pound-Drever-Hall method for stabilizing lasers to optical cavities [I3l [HI H]. 

Here we investigate two methods which utilize similar experimental arrangements, 
but which give rise to very different locking signals. The first method involves weak 
phase modulation of the probe beam, with a modulation frequency on the order of 
the natural linewidth. This method is similar to traditional lock-in detection [15]; we 
refer to this as direct modulation. In the second method, only the pump beam phase 
is modulated. Through various nonlinear mixing processes, this can generate probe 
sidebands which beat with the carrier at the detector [T6l \T7\ [T8] . This is known as 
modulation transfer spectroscopy. 

Figure [1] shows the D2 energy levels for the two most abundant isotopes of 
potassium, '^^K and ^^K, with natural abundances of 93.3% and 6.7%, respectively. The 
natural half-width at half-maximum (HWHM) is 7 = 27rx 3.017 MHz [19], which is on 
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the order of the excited-state hyperfine sphttings. Additionally, there exist so-called 
crossover resonances in saturated absorption at frequencies half-way between each pair 
of excited state transitions. These features are due to moving atoms whose Doppler 
shifts equal half the frequency difference between transitions. This leads to a group 
of six overlapping features for each hyperfine ground state of each isotope. Since the 
Doppler half-width at half-maximum, A, exceeds the ground-state hyperfine splittings 
(A = 27r X 390 MHz for ^^K at 300 K), there are also nearby pairs of ground-state 
crossover resonances for each isotope. This type of crossover resonance leads to reduced 
transmission, as the pump beam transfers greater population to the probed ground 
state via hyperfine pumping. The significant overlap between various transitions makes 
it extremely difficult to resolve individual features, and off-resonant excitation means 
cycling or closed transitions are effectively absent. In practice, the useful features are 
limited to the composite feature comprising the ^^K, F = 2 — )■ F' = 1, 2, 3 transitions 
and the one due to the two ^^K ground-state crossovers. For brevity, we will refer to 
these features as A and B, respectively (see figure [I]). These two features will be the 
focus of this work. 

Potassium is widely used in experiments with cold atoms due to the existence of 
both bosonic and fermionic species, and the availability of inexpensive laser diodes at 
the cooling wavelengths. Although we are aware of research groups using heterodyne 
spectroscopy with potassium for laser locking, to our knowledge there are no published 
examples showing what the resulting spectra should look like or how they may be 
optimized. It is therefore our aim to provide a practical study of these methods for people 
working in this field. As our interest is in laser stabilization, we attempt to optimize 
the slope of the discriminant — i.e. the rate of signal variation in V/MHz around the 
zero-crossing. This sets the gain for feedback control. In contrast to electronic gain, 
there is no inherent compromise between optical gain and total loop bandwidth, making 
it beneficial to maximize the slope optically. Despite the overlap of multiple transitions, 
we will show that narrow (< 10 MHz) dispersive features can still be observed. 

The rest of the paper is organized as follows. In section |2] we study conventional 
saturated absorption, and examine the feature amplitudes and widths as functions of 
vapour pressure. Perhaps surprisingly, the apparent half-width of saturated absorption 
features can decrease with increasing vapour pressure. This is shown to be due 
to the departure from the linear absorption regime. In section [3], we investigate 
heterodyne detection using direct modulation of the probe. The signals we obtain are 
well approximated by derivative signals. In section H] we perform modulation transfer 
spectroscopy by moving the modulation to the pump beam. This method results in 
a fiatter background and steeper slope for the ^^K F = 2 ^ F' = 3 feature, at the 
expense of eliminating the crossover. For modulation transfer, varying the pump and 
probe powers reveals information regarding the physical origins of the signal. Finally, 
in section [5l we compare the relative benefits and drawbacks of the two methods, and 
suggest some ways that they could be improved. 
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Figure 1. (Colour online.) Energy level schematics for naturally abundant potassium 
isotopes. Above: energy level diagrams for ^^K and ^^K. Hyperfine frequencies 
obtained from [20] are given in MHz, with excited states expanded by 10 x for 
clarity; F is the total (electronic plus nuclear) angular momentum. The vertical 
grey bars in the excited and ground states show the natural width ±7/(27r) and 
the Doppler width ±A/(27r), respectively. Below: transition frequencies for ^^K (red 
circles) and ^^K (blue squares). Upward-pointing features are normal transitions, and 
downward-pointing features are ground-state crossovers; excited-state crossovers have 
been omitted for clarity. The heights of normal transitions reflect the relative oscillator 
strengths. The centre frequency is taken to be the midpoint between the "^^K ground 
state crossovers. The horizontal grey bars are the natural and Doppler widths as above. 



2. Saturated Absorption 

We begin witli conventional saturated absorption spectroscopy [1]. Our experiments 
use a home-built external-cavity diode laser loosely based on the design presented in 
|21j . The diode is anti-reflection coated, with a centre wavelength of 770 nm (Eagleyard, 
EYP-RWE-0790-04000-0750-SOT01-000c|J). Measurements with a home-built fibre ring 



J Product names and part numbers are used for identification purposes only, and do not constitute an 
endorsement by the authors or their institution. 
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Figure 2. (Colour online.) Saturated absorption at varying cell temperature. From 
top to bottom, the curves correspond to temperatures of 21°C, 50° C, and 67° C. 



resonator set an upper limit of ~ 1 MHz for the laser linewidtli. The laser frequency 
is scanned by moving a diffraction grating on a kinematic mirror mount with a piezo- 
electric transducer. The elliptically elongated beam emitted from the laser is made 
approximately circular with an anamorphic prism pair before passing through a Faraday 
isolator to prevent feedback to the laser. Measurements of the power transmitted past 
a razor on a translation stage give 1/e^ intensity radii of = 1.18 ± 0.03 mm and 
Wy = 1.15±0.02 mm in the horizontal and vertical directions, respectively, at a distance 
of 0.4 m from the laser (uncertainties are from fits only). Similar measurements at 4 m 
give Wx = 1.59 ± 0.03 mm and Wy = 1.506 ± 0.006 mm. All spectra presented here were 
obtained at intermediate distances. 

Spectra were taken with a 75 mm Pyrex reference cell (Thorlabs, CP25075-K), 
containing potassium isotopes in their natural abundance. Because of the low vapour 
pressure of potassium at room temperature, we began by heating the cell and recording 
saturated absorption spectra at various temperatures. The reference cell was wrapped in 
a flexible resistive heater (Minco HK5464R14.6L12A) and a layer of thermal insulation, 
and the temperature was monitored with a sensor (Minco S665PDZ40B) pressed between 
the heater and the cell. Example spectra are shown in figure El Because the isotope shift 
is smaller than the Doppler shift, and due to the preponderance of ^^K, only a single 
Doppler-broadened background is evident in each scan. This background absorption 
varies strongly with relatively small changes in temperature. At 50°C the A and B 
features are both clearly visible. It is also possible to see the conglomeration of ^^K, 
F = 1 ^ F' transitions, which nearly overlap with the ^^K ground-state crossovers, 
causing a slightly dispersive shape near 230 MHz. 

In order to quantify the temperature dependence, we take alternating scans with 
and without the pump beam. By subtracting a purely Doppler-broadened spectrum from 
one with a saturating beam present, we obtain a single trace with a flat background 
and narrow peak (dip) for the A (B) feature. For simplicity we fit the features to 
Lorentzian functions and extract the amplitudes and widths (HWHM). The frequency 
scale throughout this work was calibrated by comparing the difference between the A 
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and B line centres with a Fabry-Perot spectrum analyser with variable mode spacing 
[22]. This yielded a frequency difference between A and B features of 224 MHz, with a 
4 MHz uncertainty estimated from variations over time and under different conditions 
{e.g., cavity length, cell temperature, pump and probe intensities). 

We model the temperature dependence as follows. For a weak [i.e., non-saturating) 
probe, the transmitted signal, normalized to the incident intensity, is given by exp(— a) 
where a is the absorption coefficient. For a two- level system, with A ^ 7, the saturating 
beam leads to a Lamb dip in the probe absorption coefficient [21 13] 



where an is the Doppler-broadened absorption coefficient, and C = 0^/(1 + is 
a Lorentzian function of detuning S from the centre of the Doppler-free feature. The 
amplitude and width of £, are functions of the pump beam intensity. We allow different 
centre frequencies for the Doppler background and C (recall figure [2]), and take as the 
value of the Doppler absorption coefficient at the resonance. Since «£) is proportional 
to the vapour pressure Pvap, we explicitly write an = oc-Pvap- The difference between 
transmission with and without the saturating beam is then 



The amplitude and width are now easily obtained, using the dependence of Pvap on 
temperature described in [23]. For the B feature, we simply take < to obtain a 
negative feature. 

The observed dependence of the feature amplitudes and widths on cell temperature 
are shown in figure [3|, together with the predictions of ([2]). The agreement between 
experiment and theory is excellent. The crossover (B) amplitude peaks first, around 
50°C, but is then compressed by the deep Doppler absorption; feature A is largest at 
62°C. Note that these temperatures depend on the length of the cell. Both features 
have similar maximum amplitudes for these data, but this may depend in general on 
the intensities of the pump and probe beams. The observed linewidths are larger than 
the natural width ~ 3 MHz of a single transition. However, due to the composite nature 
of each feature, this is not too surprising. The A feature comprises three transitions 
distributed over 30 MHz, and the B feature is due to a pair of transitions 9.4 MHz apart. 
Including power broadening, this is in reasonable agreement with our observations. Since 
our model neglects temperature-dependent broadening mechanisms such as collisions 
and finite transit times, the observed variations in width can be attributed solely to the 
exponential dependence of transmitted signal on the absorption coefficients. 

Based on the results of figure [31 we typically keep our cell around 46-48°C, 
corresponding to a maximum Doppler-broadened absorption of 35-40% for a weak probe. 
This is sufficient to produce a strong signal at each feature, with minimal experimental 
complexity. 



as = «£, (1 - C) 



(1) 
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Figure 3. (Colour online.) Saturated absorption amplitudes and widths with varying 
cell temperature. The pump power was 1.19 mW and the probe was 0.23 mW 
(both measured immediately before the cell), corresponding to 20 ± 4 mW/cm^ and 
3.9 ± 0.7 mW/cm^, averaged over the Gaussian beam profile, (a) Amplitude. Red 
circles are for the A feature (due to ■^^K, F ^ 2 ^ F' transitions), and blue triangles 
the B feature (due to the '^^K ground-state crossovers). Solid curves are fits derived 
from (b) Half- width at half-maximum. 



3. Direct Probe Modulation 

We now turn our attention to heterodyne detection, beginning with direct phase 
modulation of the probe beam. Our optical setup was similar to the one presented 
in [21] ; a simplified schematic is shown in figure HI The counter-propagating pump and 
probe beams had orthogonal linear polarizations. To obtain the phase modulation, 
we used a home-built electro-optic modulator (EOM). The LiNbOs crystal (Casix 
Optics) was Y-cut, measuring 2 x 2 x 25 mm, with anti-refiection coating on the 
input and output faces, and gold coating on the sides (Z-faces) to form electrodes. 
An operational amplifier and resonant resistor-inductor-capacitor (RLC) circuit were 
used to drive the crystal at 9.62 MHz. This frequency was chosen to be about half 
the linewidths observed in figure [3](b), as a compromise between obtaining a large slope 
and maintaining a derivative lineshape [121 [H]- The modulation was generated by a 
four-channel direct digital synthesis evaluation board (Analog Devices, AD9959/PCBZ- 
ND) which allowed adjustment of the relative phase between the local oscillator (LO) 
and the radio-frequency (RF) signal from the photodiode during demodulation. The 
transmitted probe beam was detected with an amplified silicon photodiode having 
100 kf2 transimpedance gain and 50 MHz bandwidth (Thorlabs PDA8A/M). A bias- 
T (Mini-Circiuts ZFBT-4R2GW) was used to separate the low- and high-frequency 
components of the signal. The low-frequency port was fed directly to an oscilloscope, 
allowing us to observe saturated absorption spectra. The high-frequency component 
was amplified (Mini-Circuits ZFL-500LN), mixed down with a double-balanced mixer 
(Mini-Circuits ZRPD-1), and filtered with a home-built active low-pass filter (4*^-order 
Butterworth) with cutoff frequency 200 kHz and a gain of 2. The LO power was +7 dBm. 
We found that higher LO powers increased the noise without significantly improving the 
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Figure 4. (Colour online.) Direct modulation spectroscopy. Left: simplified 
schematic. The heavy blue line shows the pump, and the thin red line the probe. 
The beams are shown separated for clarity, but overlap in the experiment. PBS: 
polarizing beam splitter; EOM: electro-optic modulator; DDS: direct digital synthesis; 
PD: photodiode; DBM: double-balanced mixer; LP: low-pass filter. Right: comparison 
of direct modulation and saturated absorption. The lower trace is the mixed-down 
signal (i.e., the output of the low-pass filter, LP), and the upper one the subtracted 
saturated absorption signal. (The upper curve has been shifted vertically for clarity.) 
The powers were 50 /iW for the probe and 250 fiW for the pump, corresponding to 
0.85 ± 0.16 mW/cm2 and 4.2 ± 0.8 mW/cm^. 

signal, while lower powers resulted in nonlinear behaviour of the phase detector at low 
RF power. 

An example of a typical signal obtained via direct probe modulation is shown 
in figure HI The signal is compared to a background-subtracted saturated absorption 
spectrum for reference. As expected, the modulation signal looks like the derivative of 
the absorption. At these lower pump and probe powers, a number of features which were 
obscured in figure [2] can now be seen. For example, the internal structure of the A and 
B transitions leads to kinks in the modulation signal. This is especially problematic for 
the crossover, as the asymmetry pushes the kink towards the zero-crossing lock point. 
We will return to this point below. As mentioned previously, the feature to the blue 
of the B feature, is due predominantly to the ^^K, F = 1 ^ F' transitions (there is a 
small contribution of opposite sign due to the ^^K ground-state crossovers). This feature 
leads to a distorted zero crossing and saturates at much lower powers than the A and 
B features. We therefore neglect it in the rest of our analysis. 

The behaviour of the direct modulation signal as a function of pump intensity is 
shown in figure [5l The slope at the zero crossing of the A feature increases until the 
saturated absorption amplitude begins to saturate and the width broadens, causing 
the derivative signal to roll over. The result is an optimum slope for the A feature 
of ~ 70 mV/MHz just below 1 mW/cm^. The B feature slope is complicated by the 
kink which occurs at low intensities, when the two crossover transitions (separated by 
9.4 MHz) are partially resolved. At 50 yuW (0.85 ±0.16 mW/cm^) the kink moves below 
zero, and then is washed out completely as the pump power increases (dashed blue 
line). The amplitude shows this transition more dramatically. The A feature amplitude 
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Figure 5. (Colour online.) Dependence of direct modulation spectra on pump 
intensity, with fixed 0.85 ± 0.16 mW/cm^ probe intensity, (a) Slope at the zero- 
crossing. The red circles (blue triangles) are data for the A (B) feature. The vertical 
dashed line shows when the kink in the B feature no longer hits the zero crossing, 
(b) Peak-to-peak amplitude, between turning points. Here and throughout this work, 
error bars reflect fit uncertainties only; when not visible, they are smaller than the plot 
markers. Not shown is an approximately 20% systematic uncertainty in intensity, due 
to variation in the Gaussian spot size away from the laser. 



saturates at much lower intensity than we expected. The saturation intensity for the 
cychng transition is 1.75 mW/cm^ [T9||20]. Assuming isotropic populations of the atomic 
Zeeman substates, this increases to 3.75 mW/cm^. Additional relaxation mechanisms 
{e.g., de-phasing or transit time) tend to further increase the saturation intensity [2]. We 
speculate that optical pumping may play a role in our observations. The importance of 
hyperfine pumping has been previously demonstrated for rubidium |25] , and the reduced 
level splittings in potassium could exacerbate this effect by increasing the of-resonant 
excitation to the non-cycling transitions. As our aim here is simply to optimize the 
locking signal, we leave the origins of this effect for future study. 

Figure |6] shows the effect of varying the probe intensity. For this choice of pump 
intensity (4.2±0.8 mW/cm^), the two features show similar slopes and amplitudes. For 
the highest probe intensities we have tested (limited by beam size and saturation of the 
photodiode), we do not see any clear saturation. 

4. Modulation Transfer 

If the phase modulation is moved from the probe to the saturating beam, we obtain 
modulation transfer spectra, shown schematically in figure [71 We found that modulation 
transfer worked better at lower frequency. An example is shown in figure [TJ with 
the modulation at 1.82 MHz. This frequency was chosen based on results in other 
labs with rubidium [211 [26]. We could not achieve a high-Q RLC circuit at this 
frequency, apparently due to the parasitic properties of available inductors. We 
therefore replaced the op-amp driving the EOM with a higher power amplifier (Amplifier 
Solutions Corporation ASC2832). We also found that the signal amplitude more than 
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Dependence of direct modulation spectra on probe 
intensity, with fixed 4.2 ± 0.8 mW/cm^ pump intensity, (a) Slope at the zero-crossing. 
Points are as in figure [5] (b) Peak-to-peak amplitude. 
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Figure 7. (Colour online.) Modulation transfer spectroscopy. Left: simplified 
schematic. QWP: quarter-wave plate; all other notations are as in figure |4l Right: 
comparison with saturated absorption. Below is the heterodyne signal, and above is 
the subtracted saturated absorption signal (shifted vertically for clarity.) The pump 
and probe powers are the same as in figure |4l 



doubled if we used a^-a^ polarizations {i.e., right-hand circular probe and left-hand 
circular pump), at the expense of nearly complete suppression of the already small B 
feature. Polarization-dependent effects have been investigated for linear pump and probe 
polarizations in ®''Rb [271 EH], but to our knowledge the a^-a^ configuration has not 
been studied. No enhancement was observed in our experiments with direct modulation. 

The modulation transfer spectrum in figure [7] differs from the direct modulation 
spectrum (figure H]) in a few important ways. As mentioned above, there is virtually no 
B feature, and the A feature is noticeably steeper, and with flatter background. The 
locking signal is well described by the derivative of a single Lorentzian, with observed 
widths (HWHM) as small as 4.51 ± 0.11 MHz. These charactersitics are all attributed 
to the nonlinear nature of the processes generating the probe sidebands. Modulation 
transfer spectroscopy is typically strongest for cycling transitions [231 [271 12H], which 
allow the inherently weak four-wave mixing of pump and probe fields to occur many 
times. This also strongly suppresses the Doppler background. 

Figure [8] shows the dependence of the modulation transfer signal on the pump 
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Figure 8. (Colour online.) Dependence of modulation transfer spectroscopy on pump 
intensity. The probe intensity was fixed at 0.85 ± 0.16 mW/cm^. (a) Slope at the zero 
crossing, (b) Peak-to-peak amplitude, for small pump intensities. 



beam intensity. The amplitude and slope grow linearly at first, eventually giving way 
to broadening as in the case of direct modulation. However, for modulation transfer the 
optimum pump intensity is much higher and the slope is an order of magnitude greater. 

The dependence of the modulation transfer spectrum on the probe intensity is 
shown in figure M The slope grows linearly and then begins to roll over. The amplitude 
is also linear for small probe intensity. Together with the linear dependence on pump 
intensity in figure |8]^a), this suggests the physical mechanism of modulated hole burning 
\16\ . The pump carrier and either sideband combine to generate a Lamb dip of oscillating 
depth. This generates a probe sideband which beats with the carrier at the detector. 
For comparison, we show the data obtained with 9.62 MHz modulation, together with 
a fit to a power law giving an exponent 2.01 ± 0.05. This behaviour is characteristic 
of reflection, which can arise in two ways. First, the pump and probe beams burn a 
standing wave in the atomic ground-state population. A pump sideband then Bragg- 
reflects in the backwards direction, creating a sideband on the probe. Alternatively, 
a pump sideband can interfere with the probe to form a moving population grating. 
In this case the pump carrier is back-reflected and Doppler-shifted, again resulting in 
a probe sideband. Third-order perturbation theory predicts that these features vanish 
under thermal averaging when A ^ 7, but higher-order interactions were shown to 
result in a signal which is proportional to pump intensity (not shown) and quadratic in 
probe intensity [I6], consistent with our observations. 

5. Discussion 

We have investigated saturated absorption spectroscopy of the D2 lines of potassium, 
and heterodyne spectroscopy using direct phase modulation of the probe and modulation 
transfer from the pump. It was shown that the saturated absorption signal strength can 
be increased by an order of magnitude or more with moderate heating of the vapour, and 
the exponential absorption of the probe beam can lead to either an increase or decrease in 
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Figure 9. (Colour online.) Dependence of modulation transfer spectroscopy on probe 
intensity. The pump intensity was fixed at 4.2 ± 0.8 mW/cm^. (a) Slope at the zero 
crossing, for 1.82 MHz modulation, (b) Peak-to-peak amplitude. Red circles are data 
with 1.82 MHz modulation, and green squares with 9.62 MHz. The solid red line is 
a linear fit for low intensities, and the green curve is a power-law fit consistent with 
quadratic growth. 

the apparent half- width of the features. The comparison between direct modulation and 
modulation transfer showed that the two methods differ significantly despite the similar 
experimental arrangements. The direct modulation signal is roughly proportional to 
the derivative of the saturated absorption, yielding frequency discriminants for both 
dominant spectroscopic features. In contrast, the modulation transfer spectrum showed 
narrower features with essentially no Doppler background. Although there was no 
crossover feature, the maximum slope near the cycling transition was significantly higher 
than could be obtained by direct modulation. Investigation of the dependence of pump 
and probe intensities suggested this feature was dominated by modulated hole burning 
at low modulation frequency, and by Bragg-reflection from atomic population gratings 
at high frequency. 

As our main interest is in laser stabilization, we should assess whether these signals 
are sufficient for our experiments. As discussed above, a larger slope increases the overall 
servo gain, which reduces the amplitude of frequency fluctuations when locked. To set 
the relevant scale, we compare our results to the tuning transfer functions of our laser: 
260 ± 13 MHz/V by piezo and 5000 ± 250 MHz/V by injection current j§| Conservatively 
taking the piezo value, this sets the overall gain (without PID loop) to be ~ 25 for direct 
modulation and ~ 150 for modulation transfer. With added gain from a PID controller, 
which can easily exceed 100 over our 200 kHz bandwidth, we therefore expect either 
method to be sufficient. We should also consider the noise level, which was typically 
~ 15 mVrms for both methods under optimum conditions. In terms of signal-to-noise 
ratio, both methods again seem sufficient. However, taking the ratio of noise to slope 
gives estimated lower bounds on the short-term frequency fluctuations of ~ 150 kHz 
and ~ 25 kHz for direct modulation and modulation transfer, respectively. Since these 

§ Our diode tunes 250±12 MHz/mA, and our current controller (Thorlabs LDC202C) has a modulation 
input with l/(50r2) transconductance. 
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estimates neglect limitations due to finite servo bandwidth and additional technical 
noise, modulation transfer may be preferred for more demanding applications. Finally, 
we consider the capture range, defined as the frequency range around the zero crossing 
over which the slope preserves its sign. This was typically about 15 MHz for direct 
modulation, and 5 MHz for modulation transfer. The optimum capture range represents 
some compromise between tightness of lock and stability against perturbations. A 
preference for large or small capture range therefore depends on the details of the 
experimental arrangement and application. 

In comparing the two methods it is clear that modulation transfer produces a 
far greater slope, but is unsuitable for locking in the vicinity of the crossover (B) 
feature. Due to the strong suppression of the Doppler background, modulation transfer 
is expected to give a more well-defined lock point, in the sense of being very near the F = 
2 — >■ F' = 3 resonance. In contrast, the direct modulation signal exhibits contributions 
from numerous transitions, and the Doppler background causes an additional shift in the 
actual lock point. However, as stated above, the small capture range of the modulation 
transfer signal, which is in many ways a virtue, may cause the laser to become more 
easily unlocked under the influence of environmental disturbances. We therefore expect 
direct modulation to be more robust under such circumstances. 

There are a number of further investigations which could be undertaken. For 
example, it may be worth iteratively optimizing the vapour temperature and pump 
and probe intensities. This could be especially useful for modulation transfer, which 
bears relatively little resemblance to saturated absorption. We were also unable to 
perform a satisfactory optimization of the modulation frequency and sideband fraction, 
which can strongly affect the signal shapes and amplitudes [T2l [24] . For this purpose, it 
would be better to have a wideband EOM and high-voltage amplifier, which were not 
available to us. One can in principle also improve the signal by expanding the beams 
with a telescope [21]. In cases when one is limited by the beam intensity in the vapour, 
but not optical power at the detector, this can increase the signal amplitude and slope. 
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